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TIME DECAY OF SCALING CRITICAL ELECTROMAGNETIC 
SCHRODINGER FLOWS 

LUCA FANELLI, VERONICA FELLI, MARCO A. FONTELOS, AND ANA PRIMO 



Abstract. We obtain a representation formula for solutions to Schrodinger 
equations with a class of homogeneous, scaling-critical electromagnetic poten- 
tials. As a consequence, we prove the sharp L 1 — ¥ L°° time decay estimate for 
the 3D-inverse square and the 2D-Aharonov-Bohm potentials. 
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1. Introduction and statements of the results 

This work is concerned with the dispersive property of the following class of 
Schrodinger equations with singular homogeneous electromagnetic potentials 
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here u = u(x,t) : R N+1 -> C, N ^ 2, a e L°°(S N - 1 ,R), § N ^ denotes the 
unit (N — l)-dimensional sphere, and A <G C 1 (S iV_1 ,M JV ) satisfies the following 
transversality condition 

(1.2) A(0)-0 = O for all 9 € S^ -1 . 

We always denote by r :— \x\, — x/\x\, so that x = r6. 

Equation (ll.lj) describes the dynamics of a (non relativistic) particle under the 
action of a fixed external electromagnetic field (E, B) given by 



E(x) := V 



1 (r)' 



B(x) = d 



' a (r) 



where B is the differential of the linear 1-form associated to the vector field 
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3, due to the identification between 1-forms and 2-forms, the 

TV = 3, 



magnetic field B is in fact determined by the vector field curl ,L\ X " , in the sense 



B{x)v = cm\(^f^) xv, 

for any vector v € R 3 , the cross staying for the usual vectorial product 
Under the transversality condition (jl.2|) . the hamiltonian 

a (r)V , «(r) 



(1.3) 
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formally acts on functions / : R w — > C as 

l A (R)| 2 +«(R)-* div S~- lA (R) A (r) 

C A ,af = -A/ + — ^ K -^— 2 ^ / - 2% —r^- ■ V/, 

M 2 M 

where divgw-i A denotes the Riemannian divergence of A on the unit sphere S^ -1 
endowed with the standard metric. 

The free Schrodingcr equation, i.e. (|l.l|l with A = and a = 0, can be somehow 
considered as the canonical example of dispersive equation. The unique solution 
u € C(R; L 2 (R N )) of the Cauchy problem 

(14) Uu t = -An 

can be explicitly written as follows: 

(1.5) u(x, t) = e ltA f(x) := -J-g-e*^ f e"^e^/(y) dy. 

This shows that, up to scalings and modulations, u is the Fourier transform of the 
initial datum /. Formula (|1.5[) contains most of the relevant informations about the 
dispersion which arises along the evolution of the Schrodinger flow. In dimension 
N = 1 , the evolution of an initial wave packet of the form 

f K (x)=e tKx e-^, KeN, 

gives an important description of the phenomenon. Inserting / = Jk hi (|l-5[> gives 
in turn the solution 

_ K 2 1 z 2 

UK\x,t) = e 2 uo(x - iK,2t); u (z,i) = e *wm . 

y/it+ 1 

This shows that each wave travels with a speed which is proportional to the fre- 
quency K, and describes both the phenomenon and the terminology of dispersion. 
The above property can be quantified in terms of a priori estimates for solutions 
to (|1.4p . A first consequence of (|1.5|) is the time decay 

M ll^/HL^^iyll/ll^. »« I + I = i, 

for some C = C(p, N) > independent on t and /. In the cases p = oo,p = 2, (jl.6[) 
can be easily obtained by (|1.5|) and Plancherel; the rest of the range 2 < p < oo 
follows by Riesz-Thorin interpolation. These inequalities play a fundamental role in 
many different fields, including scattering theory, harmonic analysis and nonlinear 
analysis. In particular, they standardly imply the following Strichartz estimates 

(1.7) ll e " A /|L?L« <C||/|| L 2, 

for some C > 0, where L\L q x :— LP (R; L q (R N ')) and the couple (p, q) satisfies the 
scaling condition 

2 N N 

(1.8) - + - = -, p^2, (p,q,7V)^(2,oo,2). 
P 9 2 

The first result in this style has been obtained by Segal in [46] for the wave equation; 
then it was generalized by Strichartz in [49. in connection with the Restriction 
Theorem by Tomas in 50 . Later, Ginibre and Velo introduced in [23] (see also 
[24] ) a different point of view which was extensively used by Yajima in [54] to prove 
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a large amount of inequalities for the linear Schrodinger equation. Finally, Keel and 
Tao in |33j completed the picture of estimates (|1.7|) , proving the difficult endpoint 
estimate p = 2, via bilinear techniques, for an abstract propagator verifying a time 
decay estimate in the spirit of (jl.61) . 

Time decay and Strichartz estimates turn out to be a fundamental tool in the 
nonlinear applications, and consequently a large literature has been devoted, in the 
last years, to obtain them in more general situations, as for example perturbations of 
the Schrodinger equation with linear lower order terms, as in (jl.ip . In particular, 
since less regular terms usually arise in the physically relevant models, a deep 
effort has been spent in order to overcome the difficulty deriving from the fact 
that the Fourier transform does not fit well with differential operators with rough 
coefficients. Among these, electromagnetic Schrodinger hamiltonians have been 
object of study in several papers. 

An electromagnetic Schrodinger equation has the form 

(1.9) iut = {-i\7 + A(x)) 2 u + V{x)u, 

where u = u(x,t) : R N+1 -> C, A : R N -> R N , and V : R N ->• R. Homogeneous 
potentials like 

(1.10) W~£p \v\~^ 

represent a threshold for the validity of estimates (|1.6|l and (|1.7|) . as shown by 
Goldberg, Vega and Visciglia in [26], when A = and later generalized by Fanelli 
and Garcia in [18] for 4^0 (actually, the authors in [TBI [2£] disprove Strichartz 
estimates, and a byproduct of this fact is the failure of the usual time decay es- 
timates). Notice that, for potentials as the ones in (|1.10|) . equation (|1.9|) remains 
invariant under the usual scaling u\(x,t) = u(x/X,t/X 2 ), A > 0, and this is why 
we refer to it as the scaling-critical situation. We also recall that equation (|1.9[) is 
gauge invariant, namely if u is a solution to (|1.9|) . then v — e % ^ x 'u solves the same 
equation, with A replaced by A + V</> and the same magnetic field B. 

In the purely electric case A = 0, several authors studied the time dispersion 
when the potential V is close to the scaling invariant case (|1.10[) (see [3] [4] [5J [T5j [22] 
[231 l2"9l I3TJ1 |4"UI US] and the references therein, both for Schrodinger and wave equa- 
tions, and also the useful survey [35] ). A typical perturbative approach consists in 
writing the action of the flow e lt< - v > via spectral theorem, and then reducing mat- 
ters of proving the desired estimate to perform a suitable analysis of the resolvent 
of —A + V, in the Agmon-Hormander style. We refer to the results by Goldberg- 
Schlag [25] and Rodnianski-Schlag [43], in which also time dependent potentials 
are treated, as standard examples of this technique for Schrodinger equations; in 
these papers, time decay estimates are obtained under integrability conditions on 
V which are close to the scaling invariant case (jl.101) , but do not include the critical 
behavior l/|x| 2 , due to the perturbation character of the strategy. Another possible 
approach consists in studying the mapping properties of the wave operators in L p , 
and obtaining the time decay for the perturbed flow e lt ^ A ~ ys> as a consequence of 
(|1.6p . via interwining properties. This point of view was introduced by Yajima in 
[55l l56l [57] and then followed by different authors (see e.g. [TTJ [52] [53]). Since it 
leads to a much stronger result, the integrability conditions which are needed for the 
potential V are usually far from being optimal in the sense of (jl.lOj) . The unique 
situation in which, at our knowledge, the L p -boundedness of the wave operators is 
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proved under almost sharp assumptions on V is the ID-case, as it has been proven 
in [11]. About Strichartz estimates for e ltlyA ~ v \ the situation is quite clear, thanks 
to the results obtained by Burq, Planchon, Stalker and Tahvildar-Zadeh in [JJ [5]; 
the authors can prove a suitable Morawetz-type estimate for the perturbed resol- 
vent, by multiplier techniques, which implies, together with its free countepart and 
free Strichartz, the Strichartz estimates for a class of potentials V which includes 
the ones which are critical in the sense of (jl.lOj) . 

The situation in the electromagnetic case A ^ 0, V ^ is quite more complicated 
and weaker results are available. Some additional difficulties, performing the above 
mentioned approach, come into play, due to the introduction of a first order term 
in the equation, which makes more complicate the analysis of the resolvent (see e.g. 
[12]). On the other hand, some results are available, both for estimates like (|1.6|) 
and (|1.7j) . under suitable conditions on the potentials in (|1.9|) . which as far as we 
know never permit to recover the critical cases as in (jl.lOp (see e.g. EH [TSJ [T7J 
[12 [33 [Ml HI1 [3S [43 E7J [48] and the references therein). 

In view of the above considerations, it should be quite interesting to produce 
a tool which might permit to prove the decay estimates (|1.6[) for equation (|1.1|) . 
in which the potentials are scaling-critical. The main goal of this manuscript is 
to give an explicit representation formula for solutions to (|1.1|) . which is in fact 
a generalization of (jl.5[) . In the approach we follow in the sequel, the critical 
homogeneities and the transversality condition (|1.2p play a fundamental role. We 
are now ready to prepare the setting of our main results. 

A key role in the representation formula we are going to derive in section [4] is 
played by the spectrum of the angular component of the operator £a,o on the unit 
(N — l)-dimensional sphere S N ~ 1 , i.e. of the operator 

(1.11) L A , a = (-iV S N-i+A) 2 + a(8) 

= -Agiv-i + (| A| 2 + a(0) - i divgN-i A) - 2i A • V s n-i . 

By classical spectral theory, La, admits a diverging sequence of real eigenvalues 
with finite multiplicity /ii(A, a) ^ /^(A, o) ^ • • ■ ^ Mfc(A, a) ^ • • • , see [321 Lemma 
A. 5]. To each fc e N, fc > 1, we associate a L 2 (§ Ar_1 ,C) -normalized eigenfunction 
ipk of the operator La, a on S^ -1 corresponding to the fc-th eigenvalue /U&(A, a), 
i.e. satisfying 

\j SN - 1 \M0)\ 2 dS(6) = l. 

In the enumeration /Ui(A, a) ^ /^(A, a) ^ ••• ^ /«fc(A, a) ^ • • • we repeat each 
eigenvalue as many times as its multiplicity; thus exactly one eigenfunction ip k 
corresponds to each index k € N. We can choose the functions ipk in such a way 
that they form an orthonormal basis of L 2 (E> N ~ 1 , C) . We also introduce the numbers 



7V-2 fN-2\ 2 ,. K „ UN-2^ 2 



(1.13) a k := — W I — 2~ ) + Mfc(A, a), fa := W I -^— I + Mfe (A, a), 

so that j3k = o 2 ~ ak ' ^ or ^ = •"■> ^' ' • ' ' which will come into play in the sequel. 
Notice that a\ ^ a-2 ^ 03 . . . . 
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N -2^ 2 



Under the condition 

(1.14) /ii(A,o)>- 

the quadratic form associated to £a,o is positive definite (see Section [2] below 
and the paper [19j). thus implying that the hamiltonian £a,o is a symmetric 
semi-bounded operator on L 2 (M. N ;C) which then admits a self-adjoint extension 
(Friedrichs extension) with the natural form domain. As a consequence, under 
assumption (11.14)) the unitary flow e ltCfiL - a is well defined on the domain of £a,o 
by Spectral Theorem; therefore, for every uq £ L 2 (R N ;C), there exists a unique 
solution u(-,t) := e ltCA °-u £ C(R;L 2 (R N )) to (JTTTJ) with u(x,0) = u Q (x). 

Remark 1.1. We notice that 



u(-,-s) := e !s£a '"u = e isC - A ^u^; 

henceforth, for the sake of simplify and without loss of generality, in the sequel we 
consider u = u(x, t) : R N x [0, +oo) -► C. 

The main theorem of the present paper provides a representation formula for 
such solution in terms of the following kernel 



(1.15) K(x,y) =^r Pk j-a h {\x\\y\)Mil\)Mfi\)> 



fe=i 



where a^, flk are defined in (|1.13|) and, for every v £ R, 

j u (r) :=r ~ J v+ e^l (r) 
with J v denoting the Bessel function of the first kind 



'.<(') =(!!£; 



{-if f^ 2h 



2) ^T{k + l)T{k + v+l)\2 

The following lemma provides uniform convergence on compacts of the queue of 
the series in (jl,15|) . 

Lemma 1.2. Let a £ L 00 ^ -1 ,!*), Ae C 1 (S N - 1 ,R N ) such that (T2) and (jlTl)) 
hold. Then there exists k<j ^ 1 such that the series 



£ i-e>j- ah (\x\\v\)Mil\)Mtit) 

fe=fc + l 

is uniformly convergent on compacts and 

K(x,y)-f2^ k J-«A\x\\y\)MTk)^^ eL ^ R2N > C y 
fc=i 

We can now state the main result of this paper. 

Theorem 1.3 (Representation formula). Leta£ L°°(S N - l ,M.), A £ C 1 (S N - 1 ,R N ) 
such that (jl.21) and (|1.14p hold. Let £a,<i as in (jl.31) and K as in (|1 . 15[) . If 

u £ L 2 (R N ) and u(x,t) = e itCA ' a u (x), then, for all t > 0, 

(Lie) *'l = ^//fes)«*** 
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Remark 1.4. The integral at the right hand side of (|1.16[) is understood in the 
sense the improper multiple integrals, i.e. 

where B R := {y e R N : \y\ < R}. 

Remark 1.5. Formula (jl.161) is in fact a generalization of (jl.5[) . Indeed, in the 
free case, i.e. A = and a = 0, the operator La. a reduces to the Laplace Beltrami 
operator — Agjv-i, whose eigenvalues are given by 

\ t = {N-2+l% £-0,1,2,..., 
having the £-th eigenvalue Xe multiplicity 

_ (N - 3 + £)\(N + 2£ - 2) 
me ~ e\(N-2)\ ' 

and whose eigenfunctions coincide with the usual spherical harmonics. For every 
£ ^ 0, let {V«, m } m =i ) 2,...,m< be a L 2 (S Ar_1 , C 7V )-orthonormal basis of the eigenspace 
of — A§jv-i associated to A^ with Y^. m being spherical harmonics of degree £. Hence 
we have that 

if fc = 1, thcn/ii(0,0) = A = 0, ai=0, ft = E f 1 , 

if k > 1 and ^„=o m » < k ^ X^=o m ™' tlien \a k = -£ , 

{/3 k = ^+t 

{4>k}k=l,2,... = {Ye 7 m}e=l,2,... 

m— 1,2... ..mi 

The Jacobi-Anger expansion for plane waves combined with the Addition Theo- 
rem for spherical harmonics (see for example [28[ formula (4.8.3), p. 116] and [6, 
Corollary 1]) yields 

oo / me 
iv — 2 f 

1=0 ^m=l 

for all x, y € ffi . Then in the free case A = 0, a = 0, we have that 

K(x,y) = sJV , w _ 2 , 

(27TJ 2 £ 2 

which, together with (|1.16|) and taking into account that if t = — s, s > 0, then 



u(-,—s) = e ls ( A )wo, gives in turn (11.51) . 

Remark 1.6. We remark that, for every y £ W N fixed, the function K(-,y) formally 
solves the equation 

£A, a K(;y) = \y\ 2 K(;y), 
as one can easily check; in fact this fits with the free case, in which K is the plane 
wave e~ tx ' v , up to constants. 

Formula (jl.161) is not present in the literature, at our knowledge; moreover, as 
far as we understand, it should provide a fundamental tool for several different 
applications. A first immediate consequence of the representation formula (|1.16|) is 
the following corollary. 
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Corollary 1.7 (Time decay). Let a € L 00 ^- 1 ^), A e C 1 (§ N - 1 ,R N ) such that 
([L2]l and ([TTTI|l /ioZd 7 and £ A ,a as m (TQ|) . If 



(1.17) sup |if(a;,|/)| < +00, 

x,y£M N 
with K as in (|1 . 15|> , then the following estimate holds 

(1.18) ||e it£A '°/(-)|| LP < ^_^ II/IIl p ', P€[2,+oo], ± + ± = 1, 

/or some C = C(A, a,p) > which does not depend on t and f . 

Proof. The proof is quite immediate. Formula (|1.16|) . ()1.17[) . and Remark [Ll] au- 
tomatically yield (jl.181) in the case p — 00. The rest of the range p ^ 2 in (|1.18[) 
then follows by interpolation with the I? conservation. □ 

Remark 1.8. Once the matters to prove a time decay estimate are reduced to the 
study of the kernel K in (|1.15p . the behavior of the spherical Bessel functions j- ak 
comes into play. The crucial fact to notice is that condition f|l .17[) is strictly related 
to the requirement oc\ ^ which in other words, by (|1 . 13|) means Hi(A, a) ^ 0. 
It is easy to verify that (jl.171) implies that /J,i(A, a) ^ 0, while arguing as in the 
proof of Lemma [L2l we can easily check that Hi(A, a) ^ implies that K is locally 
bounded. 

We are strongly motivated by the examples in the sequel to conjecture that in 
fact conditions (|1.17[) and ni(A,a) ^ are equivalent. 

We now pass to give a couple of relevant examples in which the abstract assump- 
tion (|1.17p can be checked by hands, and the optimal time decay can be obtained 
by working directly on the representation formula (|1.16|) . 

1.1. Application 1: Aharonov-Bohm field. We start with a 2D example of 
purely magnetic field, which is given by potentials associated to thin solenoids: if 
the radius of the solenoid tends to zero while the flux through it remains constant, 
then the particle is subject to a <5-type magnetic field, which is called Aharonov- 
Bohm field. A vector potential associated to the Aharonov-Bohm magnetic field in 
R 2 has the form 

(1.19) A(x 1 ,x 2 ) = a(-^,-^-\, ( Xl ,X 2 )€M. 2 , 

with a£l representing the circulation of A. around the solenoid. Notice that the 
potential in (jl.191) is singular at x = 0, homogeneous of degree —1 and satisfies the 
transversality condition (|1.2[) . 

This situation corresponds to problem (|1.1[) with 

N = 2, A(0) =A(cost, sin t) = a(- sin t, cost), a{9) = 0, 

so that equation (jl.ll) takes the form 

(1.20) ^=(-*V + «(-^,^))\, 

with x = (xi, X2) € IR 2 . In this case, an explicit calculation yields 
{ Mfc (A, 0) : k G N \ {0}} = {(a - j) 2 : 3 € Z}, 
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(see e.g. [34] for details), and in particular 

A*i(A,0) = (dist(a,Z)) 2 ^ 0. 

If dist(a, Z) t^ i, then all the eigenvalues are simple and the eigenspace associated 
to the eigenvalue (a — j) 2 is generated by tp(cost, sini) = e -2 - 7 *. If dist(a,Z) = i, 
then all the eigenvalues have multiplicity 2. The following result is an interesting 
consequence of Theorem 11.31 

Theorem 1.9 (Time decay for Aharonov-Bohm) . Let N = 2, a € M., and define 

^■■={-^ + a[- WW 

Then the following estimate holds 

(1.21) \\e itC "f(-)\\ LP < J! ry ll/ILy, pg[2,+oo], 1 + ^7 = 1- 

\t\ y? p) P P 

for some constant C — C(a,p) > which does not depend on t and f. 

Remark 1.10. Since curl^4(x) = if x ^ 0, the action of the magnetic field 
in the Aharonov-Bohm case is concentrated at the origin. However the poten- 
tial A cannot be eliminated by gauge transformations; this is in fact the peculiar 
property of Aharonov-Bohm fields, which indeed describe an interesting difference 
between the classical and quantum version of the electromagnetic theory. Due to 
the above remark, estimates (|1.21[) are not trivial, as far as we understand, and at 
our knowledge they are not known. We finally stress that the algebraic structure 
of Aharonov-Bohm potentials is exactly the one which has been used in [T5] in 
order to disprove the dispersion (in that case Strichartz inequalities) in the case of 
magnetic field which decay less than |x| _1 at infinity, in dimension N ^ 3. In 2D, 
counterexamples as the ones in |18j are missing, and it is still unclear what might 
happen for potentials with less decay than the one in (|1.20|) . 

1.2. Application 2: The inverse square potential. We also present an ap- 
plication of formula (|1.16[) in the case of perturbation of the Laplace operator in 
dimension N = 3 with an inverse square electric potential; more precisely, we con- 
sider problem (|1.1|) with A = and a(9) = a = constant, so that the hamiltonian 
(|1.3|) takes the form 

(1.22) C a :=-A+-^, inM 3 , 

\x\ 

and condition (11.14)) reads as a > — i. Since in this case the angular eigenvalue 
problem (jl.lip becomes 

q 23} / -A§2ipk = (pk(0,a)-a)ip k , in S 2 , 

we have that {'/'fc}^! are the well known spherical harmonics and 

yK fe (0,a) = a + /i/c(O,0) . 

As in Remark [HI for every i ^ 0, let {Y^ m } m=ia ..... 2 i+i be a L 2 (S N - 1 ,C N )- 
ortho normal basis for the eigenspace of — A§2 associated to the ^-th eigenvalue 
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A^ = £(£ + 1) of — A§2 (which has multiplicity me — 2£ + 1), with Ye <m being 
spherical harmonics of degree £. Hence we have that 



(1.24) if k = 1, then /Ui(0, a) — a, oi\ = \ — J\ + a, 

(WO, a) = a + £(£+l) = 

(1.25) if k > 1 and £ 2 < k < {£ + l) 2 , then <^ / j . 

W = \-^(£+\) +a 

Notice that a,\ ^ if and only if a ^ 0. We define the well known zonal functions 

(1.26) Zf\e') = J2 Yi, m {0)Y t , m (e'), 9,9' e^ 2 , 1 = 0,1,2,.... 

m— 1 

The study of the kernel in (|1.15p , which can be rewritten as 

oo 

(1.27) ^ y )=E^^ TT7 ^^, +V (iTT7^(l-IH)^(§)' 

1=0 

permits to prove the following result. 

Theorem 1.11 (Time decay for inverse square potentials). Let N = 3, a > — h, 

and define L a by (|1.22p . 

i) If a ^ 0, then the following estimates hold 

(1.28) \\e ltCa f(-)\\ LP ^ ,,? o II/IIl*'. Pe[2,+oo], - + ^ = 1, 

\t\ 3 (2-p) V P 

for some constant C = C(a,p) > which does not depend on t and f. 
ii) If — i < a < 0, let ct\ as in \1.2J$ , and define 

\\u\\ p , ai :=(j 3 (l + \x\-^) 2 - p \u(x)\ p dx\ , p>l. 
Then the following estimates hold 

(1.29) ||e rt£ V(0|L ai < C(1+ Jf 1 i r 7 |l/llp-,a 1 , P>2, ^ + ^ = 1, 

for some constant C — C(a,p) > which does not depend on t and f. 

Remark 1.12. As far as we know, the best dispersive results concerning this 
kind of operators are about Strichartz estimates, and have been obtained by Burq, 
Planchon, Stalker and Tahvildar-Zadeh in [318]. As a fact, estimates (|1.28[) imply 
the ones obtained in [7], by the standard Ginibre-Velo and Keel-Tao techniques in 
[23l [33] . On the other hand, in [8] the authors can treat more general potentials 
with critical decay, including e.g. the cases in which a — a(x/\x\) is a 0-degree 
homogeneous function; in addition, we think that the restriction N = 3 in Theorem 
11.111 is not in fact a relevant obstruction. We are motivated to claim that a deeper 
analysis of formula (jl.161) should permit to prove the analog to Theorem ll.lll in the 
more general case a = a(x/\x\), but this will not be treated in the present paper. 

Moreover, notice that ct\ > 0, in the range —\ < a < 0, see (11.241) . so that the 
decay in (|1.29p is weaker than the usual one. We find it an interesting phenomenon, 
since on the other hand the usual Strichartz estimates are still true in this range, as 
proved in [7J. Estimates (|1.29l) are presumably sharp and, at our knowledge, new. 
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The rest of the paper is organized as follows. In Section [2j we describe the 
functional setting in which we work, in order to prepare the proof of the main result, 
Theorem 11.31 Section [3] is then devoted to the study of the spectral properties of 
a magnetic harmonic oscillator with inverse square potential, denoted by Ta, (see 
formula (|2.1ip ). which comes into play when a suitable ansatz (formula (|2.9|) ^ is 
stated; finally, Section |4] is devoted to the proof of Theorem II. 3[ while in the last 
Sections [5] and |6] we prove the applications, Theorems II .91 and II. Ill 

2. Functional setting 

Let us define the following Hilbert spaces: 

• the space T-L as the completion of C^° (R N \ {0}, C) with respect to the norm 

H\\n = ( j^ (W(x)| 2 + (M 2 + pj) |<Mx)| 2 ) dx^ 

• the space H as the completion of C£°(R , C) with respect to the norm 



U\\h={ I (|V#r)| 2 + (M 2 + l)|#r)| 

V Jon \ 



1/2 

2 Ux' 



• the space W.A as the completion of C£°(M. N \ {0},C) with respect to the 
norm 

,1/2 



R N 



(\V A <fi(x)\ 2 + (\x\ 2 + l)\cf>(x)\ 2 ) dx 



with V A </> = V0 + % A{x /\ xl U - 

\x\ 

From the above definition, it follows immediately that 

(2.1) T-L <^-> H with continuous embedding. 

A further comparison between the above defined spaces can be derived from the 
well known diamagnetic inequality (see e.g. 36 ) 



(2.2) |V|0|(aO| < 



_ , . . .A(x/\x\) ,, . 



N>2, 



which holds for a.e. xel" and for all <f> G C£°(R N \ {0},C), and the classical 
Hardy inequality (see e.g. [20 ] [27 ] ) 

\ 2 r \M~\\2 



,2,^/^-2Y f \<Kx)\< 



(2.3) / \Vc/>(x)\Ux> — r — / lz ^dx, 
Jr« V 2 / J& N \ x \ 

which holds for all <j> G C%°(^L N ,C) and N ^ 3. We notice that the presence of a 
vector potential satisfying a suitable non-degeneracy condition allows to recover a 
Hardy inequality even for N = 2. Indeed, if N = 2, (11.21) holds, and 

1 f 2 " 

(2.4) $ A := — / a(t)dt$Z, where a(t) := A(cosi, sin t) • (— sin t, cos £), 

2tt Jo 

then functions in C£°(1l N \ {0},C) satisfy the following Hardy inequality 

(2.5) (min|fc-<I>A|) 2 / ^f- dx < / 
being (minfegz |& — $a|) the best constant, as proved in 



Vu(l) + i ^ — j L lt(x) 



dx 
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Combining (|273]> . fj272]) . and ([275]) . it is easy to verify that if N > 3, then i? = 
"H = Ha, being the norms || • \\h, \\ • \\n an d || • \\-h a equivalent. If N = 2 then 
UC\H;on the other hand, if TV = 2 and ([112]), ([23]) hold, from ([22]) and ([23]) we 
deduce that % = "Ha, being the norms || • \\-h, || ■ \\-h a equivalent. 

From (|2.ip and [32] Proposition 6.1], we also deduce that 

(2.6) % is compactly embedded into L P (R N ) 

for all 

'2* = jf^2, if TV ^3, 

^+oo, if JV = 2. 

The quadratic form Qa.o associated to /3a, a, i-e. 

(2.7) Q A , a :^ 1 ' 2 (K Ar ,C)^M, 

a(x/\x\) 



2<p< 



QA,a(^) 



Va0O*O 



with PPQR^, 
norm 



C) being the completion of C° 



\x\' 



-\mv 



dx, 



\ {0},C) with respect to the 



\u\l 



— , ,i2 \u(x)\ 2 . , 

Vu(x)\ + ' , ^ ) da; 



1/2 



is positive definite if and only if (|1.14[) holds, see [T9l Lemma 2.2]. In particular, 
assumption (11.14[) ensures that the operator £\,a is semibounded from below, self- 
adjoint on L 2 with the natural form domain, and that there exists some constant 
C(N, A, a) > such that 



(2.8) 



,2 o(x/|x|) ,, ,„ \X 



\<t>W + ^\<Kx)\* 



dx>C(N, A, a) \\4>\\ 2 n , 



for all</>e"H (see [19]). 

Up to a pseudo conformal change of variable, see |32j . equation (11.11) can be 
rewritten in terms of a quantum harmonic oscillator with the singular electromag- 
netic potential, as stated in the following lemma. 

Lemma 2.1. Let [tt® hold and u € C '(R; L 2 (R N ')) 6e a solution to fTTT]) . 77ien 
(2.9) 



(?(»,*) = (l + < 2 )T" M ( v /T+f2 x ,i) f 



satisfies 

and 
(2.10) 



¥ j€C(R;X 2 (R jv )), p(ar,0)=u(a;,0), 

II^(-^)IIl 2 (r«) = ll""(-,t)|U2(RW) /or a/H € R, 



if(M) = 



(1 2 s j^a.^O^ + tHVOM) 



A representation formula for solutions u to (|1.1[) can be found by expanding the 
transformed solution ys to (|2.10[) in Fourier series with respect to an orthonormal 
basis of L 2 (R N ) consisting of eigenfunctions of the following quantum harmonic 
oscillator operator perturbed with singular homogeneous electromagnetic potentials 



(2.11) 



t a ,« :«->«*, T A , a = c A . a + -\x\ 2 
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acting as 

(2.12) n*(T A , a v,w)n 



a(# 

V \x 



|.,| 



Va«(i) • Vaw(x) — — v(x)w(x) H — v(x)w(x) dx, 

\x\ z 4 J 

for all v, w £ H, where W denotes the dual space of % and u* (', -)u is the corre- 
sponding duality product. 

3. The spectrum of T A , a 

From (I2.6[) . (|2.8p . and classical spectral theory, we can easily deduce the following 
abstract description of the spectrum of TA, a - 

Lemma 3.1. Let A £ C 1 (S N - 1 ,'R N ) and a £ L 00 ^- 1 ) such that (PHl holds. 
Then the spectrum of the operator T A , a defined in h2.11\WAty) consists of a diverg- 
ing sequence of real eigenvalues with finite multiplicity. Moreover, there exists an 
orthonormal basis of L 2 (R N ) whose elements belong to H and are eigenf unctions 

0fTx,a- 

The following proposition gives a complete description of the spectrum of the 
operator Ta,o- 

Proposition 3.2. The set of the eigenvalues of the operator Tj±. a is 

J7 m:fc : k,m£ N,k ^ l} 

where 



N N-2 //iV-2 N2 



(3.1) 7 m , k =2m-a t + -, a k = — \ I — - — J + ^ fe (A, a), 

and /Xfe(A, a) is the k-th eigenvalue of the operator Lj±,a on the sphere S^ -1 . Each 
eigenvalue J m ,k has finite multiplicity equal to 

and a basis of the corresponding eigenspace is 

N 

~2 
where 

(3-2) V n ,j(x) = | JC |-«ie- J ^P i ,„(J^)^(^), 

ipj is an eigenfunction of the operator La, a on the sphere E> N ~ X associated to the 
j-th eigenvalue /ij(A, a) as in (|1.12p . and Pj >n is the polynomial of degree n given 
by 

(-n)i t l 



V n .j :j,n£ N,j > 1, 7 m , fe = 2n - a } ■ + — }■ , 



Pj,n(t)=J2 



o(f-^)^ r 



ji-1 , 



denoting as (s) i; for all s £ R, i/ie Pochhammer's symbol (s).j = Jli=o( s + •?')> 
(s)o = l- 
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Proof. Assume that 7 is an eigenvalue of Ta_,o, and g € H\ {0} is a corresponding 
eigenfunction, so that 

(3-3) f-iV+^Sl) g{x) + a JM g{x) + K g{x)=jg{x) 

in a weak "H-sense. From classical elliptic regularity theory, g G C l( ^"(M. N \ {0}, C). 
Hence g can be expanded as 

00 
g(x) = g(r$) = £ fc (r)ifo(0) in i 2 ^" 1 ), 

where r = |x| € (0, +00), = cc/ j cc j G S^ -1 , and 



0fc(r)= / g{rd)^ k {9)dS{9). 
Equations (|1 . 12|) and (J3T3J) imply that, for every k, 

ui , N ~ 1 JJ , ( Vk 

<Pk + 7 - — , 

r \ r z 4 



(3.4) ^ + __^ fc+ ( 7 _^__)0 fc = o in(0,+oo). 



Since g £ "H, we have that 



(3.5) 00 >/ g 2 (x)dx = f ( f g 2 {rB) alS^))^- 1 air 

Jm N Jo \J§n-i 



i) 



and 



(3.6) 00 >/ 9 -^-dx> [°° r N ^<&(r)dr. 

Jr n \ x \ Jo 

For all k = 1,2, . . . and £ > 0, we define «;&(£) = (2£)~2~e20fc(v2l), with a*, as in 
(|3~Tj) . From (|X4| . u> fe satisfies 

twUt) + (f - « fc - t) w' k (t) -(^-^-T\ Wk ( t ) = Q in (0, +00). 

Therefore, w k is a solution of the well known Kummer Confluent Hypergeometric 
Equation (see [5] and [37] )■ Then there exist A k ,B k £~R such that 

/TV cut- 'Y ./V \ /iV Qit- t iV \ 

Wk (t)=A k M(--^-±,--a k ,t)+B k u( K --^--j,--a k ,t) 1 te(0,+oo). 

Here M(c,b,t) and, respectively, U(c,b,t) denote the Kummer function (or conflu- 
ent hypergeometric function) and, respectively, the Tricomi function (or confluent 
hypergeometric function of the second kind); M(c, b, t) and U(c, b, t) are two linearly 
independent solutions to the Kummer Confluent Hypergeometric Equation 

tw"(t) + (b- t)w'{t) - cw(t) = 0, t£ (0, +00). 

Since (^ — a k) > E from the well-known asymptotics of U at (see e.g. [1]), we 
have that 

U( — - -, a k ,t] ~ const t 1 ~^ +ak as£^0+, 

I 4 2 2 2 J 



oo 
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for some const ^ depending only on iV, 7, and ak- On the other hand, M is the 
sum of the series 

' (c) n t" 

(b) n n! 

We notice that M has a finite limit at + , while its behavior at 00 is singular and 
depends on the value -c = -f + !f + § . If -f + ^ + 1 = m € N = {0, 1, 2, • • •}, 
then M(^ — ^ — ^, ^ — afe, i) is a polynomial of degree m in i, which we will 
denote as P k m , i.e., 



n=0 



P fe , ro (t) = M(-m,f -a k ,t 



„-0 V 2 afc Jn 



^(K-~.\ n |- 



If ( — t + if + 2 ) ^ ^' * nen fr° m the well-known asymptotics of M at 00 (see e.g. 
[5]) we have that 

, r { N at- 7 AT \ f jv , °fc 2 

My- -, — - a k ,t J — const e t 4+2 2 as t -> +00, 

for some const ^ depending only on AT, 7, and a/c. 

Now, let us fix k € N, fc ^ 1. From the above description, we have that 

w k (t) ~ const B k t 1 ~% +ak as t -> 0+, 

for some const 7^ 0, and hence 

2 

<^fe(r) = r _Q! *e _ T w fe (— ) - const B fe r" (JV " 2)+Qfc as r -> 0+, 

for some const 7^ 0. Therefore, condition (|3.6[) can be satisfied only for B^ = 0. If 
^ + I - f g N, then 

uife(t) ~ const A fc e*t"~ + "2 2 as i -> +00, 

for some const ^ 0, and hence 

2 
^M == r~ ak e~~Wk( — ) ~ const j4fcr~~~ 7 e r ' 4 as r — > +00, 

for some const 7^ 0. Therefore, condition (|3.5[) can be satisfied only for A k = 0. If 

2 2 

^ + ^-^j=m<GN, then r _Qfc e _ ^Pfc,m(^-) solves (|3.4[) ; moreover the function 

^.fcCs) - M~ Qfce ~^,™(^)v^(A) 

belongs to "H, thus providing an eigenfunction of L. □ 

Remark 3.3. If a{0) = 0, A = 0, the spectrum of Lo.o is described in Remark 1 1.51 

JV 

2 

tor -A + J 4 

Remark 3.4. It is easy to verify that 

if (mi,£;i) 7^ (ni2) fe) then V^,^ and Kn 2 ,fe2 are orthogonal in L (1 ) 
By Lemma l3~Tl it follows that 



so that the spectrum of T^o is ^ +N. Hence, in this case we recover the eigenvalues 

I i2 

of the Harmonic oscillator operator —A + J-J- (see e.g. [32]). 



is an orthonormal basis of L 



1 2fm,N\ 
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Remark 3.5. Denoting by I^(t) the generalized Laguerre polynomials 



* — ' \m — n I 

>i— n \ / 



m + a\ t n 
I' 



and Pk = V (^; ) + Mfc(A, a) so that 7 m ,fc = 2m + /3/- + 1, we can write 

From the well known orthogonality relation 

x a e- x LZ(x)L^(X)dx = r (» + a + l) 
o n\ 

where S n>m denotes the Kronecker delta, it is easy to check that 

4. Proof of Theorem 11.31 

The proof of Theorem 11.31 uses Lemma 11.21 which is proved below. 
Proof of Lemma 11.21 From Theorem IA.1I and Lemma IA.2I in the Appendix, 
we deduce that there exist some k Q £ N, and C, > 0, i — 1,2, such that for every 
k > ho, 

(4.1) -a k >C 1 k 5 \ 
and 

(4.2) 1^(0)1 < C 2 k &2 for all 6 £ S* -1 , 

with <5i = j±^ and 5 2 = j^zil^f 1 }- From ffD and 63> it; follows that, for 
all k > k , the fc-th term of the series K k (x,y) = i~& k j-a h (\x\\y\)ipk(^)ipk(^) 
belongs to L^ C (M. 2N , C). Furthermore, if we fix some compact set K (s K , there 
exists fcg such that | ._ ^L | < \ for every x,y £ K, and k > k' . Therefore, for 
all x, y £ K, and k > k' , we have that 

|x fc (x, y )| ^ \j- ah (M\y\)\ Mjt\)\\Mw\) 

" V 2 J A, o r(m + l)r(m-a fc + ^ + l) V 2 



r(-a fc + f) V2(-a fc + f) 



< C"fc 2& [ - I = M fe 



2 



where C" depends on £ but not on k. Weierstrass M-test and convergence of 
Ylk M k yields then the desired uniform convergence. D 

We are now ready to prove our main result, the representation formula given by 
Theorem O 
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Proof of Theorem 11.31 Let us expand the initial datum uo = u(-,0) — </?(•, 0) 
in Fourier series with respect to the orthonormal basis of L 2 (M. N ) introduced in 
Remark 13.41 as 



(4.3) u = V" Cm,kV m ,k vn. L 2 (R N ) , where c m ,k= / u (x)V m:k (x)dx, 



m.feeN 



and, for t > 0, the function <p(-,t) defined in (|2.9|) as 
(4.4) tp(;t)= E Vm,k{t)V™,k inL 2 (l 



,N\ 



m,k& 
k>l 



where 



<Pm,k(t) = / ¥>(s>*)K»,*(aOda:. 

Since (fi(z,t) satisfies (J2.10J) . we obtain that <p m ,k G C 1 (M,C) and 

*<Pm,fc(*) = ]_ "'^ 2 ym,fcft) ; Vm,fe(0) = C m: fc, 

which by integration yields <p m ,k(t) — c m ^e~ t ' rm,k arctan *. Hence expansion (|4.4|l 
can be rewritten as 

<p(z,t)= J2 c m , fe e-^^ arctant y m ^(z) inL 2 ^), for all* > 0. 

m,fcGN 
fe>l 

In view of (|4.3j) . the above series can be written as 

<p{z,t) = J2 e-^.* arctan *( / «o(y)V m , fc (»)dy , )y mifc («), 

m,fc€N VJK™ / 

fc^l 

in the sense that, for all t > 0, the above series converges in L 2 (R A '). Since tto(y) 
can be expanded as 

oo 

My) = M\y\ ft) = E ^(lyl^-Cft) m l 2 ^- 1 ), 

3=1 



where Uo,i(|2/|) = Jsjv-i u o(l2/|0)V , j(0) dS(6), we conclude that 



p — ilm k arctan £ / /*oo 2 

(*>*) = E iitT ||2 V m , k (z)( / «o, fe (r)r Ar - 1 -^P fe , m (4)e-^dr 

V 771 / w 






E 



00 g-i(/3 fc + l) arctan t 

E^(^) 2 ^r(i + /3 fe ) 



cilm arctan £ 
m— 



^ 2 _L I ~ I 2 







rz 



By [5] we know that 

„2, rvi i a. A 2 „. /■<» 



My) ■ r$S^*-**f '-*~*-'*<^* 
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where Jp k is the Bessel function of the first kind of order (3k- Therefore, 



ip(z,t) = 4^Vfc( 



fc=i 



z) r(i + /3 fc ) 

l z l ' g*(^fc+l) arctant 



E 



m + /?fc\ e 



uo,fc(r) 



. m— 

_, / /-OO /-0O 



-i2m arctan £ 



e * 



m 7 (T(l + P k + to)) 2 ' 

-s 2 -8' 2 f o J\2m+0 h +l v 



k |r«| Q *+^ vo Jo 

x J^ {y/2rs)Jp k (V2|z|s') ds ds')^ -1 ^ 



4 E^(ifr 



? — i(/3fe + l) arctan t 



x e 



fc=i 

i (arctan £+-j)/3fc 



{ss'f 



uo^lrzl-^-^e^^-^^x 



cS 2 + S /2 



_°°. ( i Vm — i(arctani+-^-)(2m+/3fc) 

E rn 7 .,, ■ I , ^ (^) 2m+fa 



rn— 



"(1 + m)r(l + /3 fc + m) 



Jp k (V2rs)JpjV2\z\s')dsdAdr 



Then, since £ (-l)m elgg^ggj ( s /)^ = j^^-^rctant+f),^ 

m— 

have 



(4.5) p(z,t) 

CX) 

= 4 E^: 



= i(/5fc ^-—arctan i) 



fc=l 



jV-2 r"° + |* 



Uo,k(r)|r2| 2 e * X fe)t (r, |z|)r dr 



where 



Z M (r,|z|)= / / SS 'e- s2 - s ' 2 J 0k (2e- l ^ ct ™ t+ ^ SS ')J Pk (^r S )J h {^\z\s')dsds'. 



From [53 formula (1), p. 395] (with t = s',p = l,a= \/2\z\, b = 2e-' l ^ Ictant+ ^s 1 
v = /3k which satisfy §?(i/) > —1 and | argp| < j), we know that 

s'e- s ' 2 J /3fc (2e- j(arctant+ 5) ss ')J /3fc (V2\z\s') ds' 



1 | 2 |2 +2c -i(2arctant + x) s 2 

r 2 - 



/9* I _i(a 



V2|z|s 



« V t*rctan£+-^) /' 

where Ip k denotes the modified Bessel function of order j3k ■ Hence 
Zk,t( r >\ z \) 

1 f°° 2 , z ,2, 2e -i(2 arctan t + 7r) s 2 

se~ s Jp k (V2rs)e a I^(V2e" l(arctant+ ^ ) |z|s)ds 



2 Jo 



1 



83 ^ (rs)e- 



|2 1 -i(2 arctan t + vr) 2 



Ip k (e~ i(arctant+ ^\z\s)ds. 
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Since I v (x) — e~^ vm J v {xei l ) (see e.g. [3] 9.6.3, p. 375]), we obtain 

l k , t (r,\z\) = Y^ l ^l S e-^^ ( — + ^\j, k (rs)J Pk {e—^ ^\z\s)ds. 

Applying [5TJ formula (1), p. 395] (with t — s, p 2 = i+s ^ , a = r, 

b = e -*arctant| z | ; v = fa w hi c h satisfy Ht(u) > -1 and | argp| < f ) and [2, 9.6.3, 
p. 375], we obtain 

(4.6) I k ,t{rM) 

1 g, i.|2 1 r 2 + | z | 2 e - 2iarctant / r | 7 |„-iarctant 

= i e -^^ e --T- e 2 (i + e-*(2«<^ '+»')/. | LilLr 

A 1 I p— z(2 arctan i+7r) \ 1 + p— i(2 arctani+7r 

1 ,.,a 1 r 2 + M 2 e -2iarctant / „|„| p --tarctant 

_ f e -/3fe« e -J4- i e - 2(1+e - i( 2a,c t .„ t+I ) ) j (j nfig 

/I 1 I p— z(2 arctan i+7r) Pfc I i _i_ _ i(2arctan t+7r) 

Noticing that 

-iarctani _ l V' + Z ) 1 _ ^ ~ l 

y/l _|_ £2 ' 1 -|- e -i(2arctant+7r) — 2£ ' 

from (I4.5[) and (|4.6[) we deduce 

(4.7) p(«,t) 

_ — zarctant °° /-co / \ 2 / . , \ 

= ! Vt^JAle-^f / ii »i r ) c V(l- , l -.. 2a L,HJ , 

J +e -i(2arctant+7r) Z^ ^ fe V|z|/ e ' ' ' 

k— 1 



JV— 2 

o Irzl 2 



M 2 (i i 2e -2iarctant N / „ I _ I .-i arctan t \ 



^ fc I X + g—«(2 arctan t+7r) I 



^^EMrK* 



2ti 

fe=i 



wo,fc(r) _x?_ _^P T fr\z\y/l + t 2 > Ar _ [ , 



/o |rz| 2 
From (|4~7) and ([2~9|l we get that for i > 0, 

(4.8) u{x,t) = (l + t 2 )~%e'*f^ip(-fJL= s ,t 



'' l 2i 






fc— i 



-/,' 



Notice that, by replacing J by J in (|4.8|l one obtains the series representation 
of the solution u R (x,t) with initial data Uq,r(x) = X-R( a; ) u o(a;) with X-r( x ) the 
characteristic function of the ball of radius R centered at the origin. Since the 
evolution by Schrodinger equation is an isometry in L 2 , we have that for all t G K 

\\u - u R \\ L 2 (t) = \\u - u 0tR \\ L2 ->• 0, as R -> oo. Hence u(-,t) = lim^^oo u R (-, t) 
inX 2 (R N ). Since 



«o,*(r)= / Mr0)M0)dS{0), 
and, by hypothesis, the queue of the series 

fe=i ^ 
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is uniformly convergent on compacts, we can exchange integral and sum and write 

e iti y_2 j / _v r 2 

u R {x,t)= \x\ 2 / / u {r9)r^e «< x 



fc=i ^ 

i|a;| 2 



dr dS(0) 



/«fe *>'*"»<»>*■ 



Letting i? — >• oo, we obtain ()1.16j) thus completing the proof of Theorem 11.31 D 

5. Proof of Theorem 11.91 

In view of Remark ll.il it is enough to prove the stated estimate for t > 0. Moreover, 
thanks to Corollary 1 1.7[ it is sufficient to prove condition (|1.17j) . namely, uniform 
boundedness of 

7v/ \ X ^ — i\a — i\— -— ii arctan — ii arctan — T /i 11 i\ 

K{x, y) = ? — rj 2^ e *'" Jl 2 e 'ie 3 « J|«-j| (l^llz/l) 

which can be written as K(x,y) = , 2 \^ Vt / ( arctan ^ — arctan^ 2 -, |x||y|) where 

W(z,s) = J2e- ila ~ Me ~ iJSJ \<*-J\(z)- 

jez 

Notice that 

, ., \oi-j, ifj<a, 

\ a ~ J\ = \ 

[J -a, if j > [a] = j , 

so that we can write 

W{z, s) = r a J2 i j e~ ijs Ja-j (z) 

j<-bo|+i 

bo |+i _ 

+ J2 e-^-^fe^ s J| Q _ J |(z) + z Q J2 i~ J e- ija Jj- a (z) 
-bo|+i i>bo|+i 

= r a 5i(z, s) + S 2 {z, s) + i a S z {z, s). 

S2(z,s) is clearly bounded. By using identity 9.1.27 in [2], 

J'v( z ) = 2( J "-i( z ) - J ^+i( z ))' 
we can compute 

(5.1) ±S 1 (z,s) = l Yl i- j e ijs (Ja +j -i(z)-J a+j+ i(z)) 

j>bo|+i 

= l^ r (, + i) e .(^ 4+j(z) _l Y i-^e^-^J a+j (z) 

J>\jo\ J>b'o|+2 

= l(i- 1 e is J aMjo \ + i(z) + J aMjol+2 (z))r^e i ^ s -(ico8s) &(«,*), 
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(5.2) ±S z {z,s) = \ J2 i- j e- ijs (J- a+j -i(z)-J- a+j+ i(z)) 

i>liol+i 

3>\jo\ J>lio|+2 

= \ (.-^-^J-o+i^i+i^) + J- a+ \j \+2(z)) i-We-tW' - (i coss) S 3 (z, s). 

Defining 

F(z, s) = i- a S 1 (z, s) + i a S 3 (z, s), 

we deduce from (|5.ip and (J5.2I) that, for every s, F(-,s) satisfies the differential 
equation 

(5.3) —F(z, s) + (i cos s) F(z, s) = g(z, s), 
with 

g(z,s) = ^- (i- l e is J a+lJo]+1 (z) + J a+M+2 (z)) r^-^ e K\M+^ 



+ y (i- l e- is J_ a+lM+1 (z) + J_ a+lj0l+2 (z)) r'*'l- 1 e- < (l«l +1 >. 
By integration of (|5.3[) we obtain that 



F(z, s) = e~ ,zcoss I F(0, a) + / e" coss <?(z', s)dz' 

Since |jo| + 1 ± a > 0, by the asymptotic behavior of Bcssel functions close to the 
origin (sec formula 9.1.7 in [2]) 



1 fx\ v 
we conclude that F(0, s) = 0, and hence 



■ M *> = nTTT) (f ) 



F(*,a) = / e-^- z '> coss .g(2',s)dz'. 

Uniform (in s and z) boundedness of F(z, s) follows from uniform boundedness of 
the function 

f(z,s) = j Z e lz ' coss {i- x e ls J a {z') + J a+l {z')) dz' 
Jo 

for any a > 0. In order to prove it, we use the identity 

2 / <77T 7T S 



(5.4) j ff ( z ) = ^_ C0S ^_^[_E) + f .( z ) 

with 



(5-5) |&(*)| < - 



2:2(1 + 2)' 
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which is a simple consequence of the asymptotic behavior of Bessel functions at 
infinity (see formula 9.2.1 in [2]). Therefore, 



/(*,*) = 



l'^"'-(^"^'-T-l 



sm z — 



— - ))dz' 



Az cos s ( ■ — 1 As 



(i- l e™&(z') + & + i(z'))dz' = h(z) + I 2 (z). 



By (|5.5p . I-2{z) is uniformly bounded. We notice now that 

& <.'«». ( e -< f e i S cos ( Z ,_™_?L\ + sin (V _ £E _ 1 

\l ttz' V V 2 47 V 2 4 

+ l) e l ( z '-^-T) + ( e « _ 1) e -<(*'-¥-i)^ 



/ is' COS S / / 26' 

27V ^ e (!( 



and since 



I 



Z Az' (COS 3 + 1) 



2 p i2'(coss- 1) 



and 



dz 



dz 



1 



1 



Vcos s + 1 7o 



2(coSS+l) giy 



e 



J /^(l-coss) p -jj/ 



Vcos s + 1 



<c, 



VI — coss Jo 

V s - 1 



Vv 



dy 
dy 



VI — COS; 



<C, 



we conclude that 7i (z) is uniformly bounded. 

Therefore, K{x 1 y) is uniformly bounded and then inequality (|1.21|) follows by 
Corollary O □ 

6. Proof of Theorem 11.111 

In view of Remark 11.11 it is sufficient to consider the case t > 0. Let TV = 3, 
a > — j , and K as in (|1.27[) . The proof of the theorem will follow from the following 
estimates for K: 



3.1 



if a.\ < 0, then sup \K(x, y)\ < +oo, 

x,yeK 3 

if a\ > (J, then sup j. — .. .. — 

x,yeR3 1+ (\x\\y\j Ul 



< +00. 



(6.2) 



Before proving the above estimates, let us show how (|6.ip and (|6.2[) imply estimates 

(jl.281) and (|1.29[) respectively, thus proving Theorem II. Ill 

We notice that if ot\ =0 then a — and there is nothing to prove since in 

this case the result reduces to classical decay estimates for the free Schrodinger 

equation. 

If o.\ < 0, in view of (|6.ip estimate (|1.28[) directly follows from Corollary [TT7J 
If a\ > 0, then (|6.2[) and (|1.16|) imply that, for some C > (independent of x 

and t), 



C 
\u(x,t)\ < — 

t2 



1 



- ai \y\- ai 
c 



\uo(y)\dy 



c .. .. , c i r \My)\, 
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for a.e. x € K 3 and all t ^ 0, which implies 



(6.3) 



\x\ ai l+t ai 

+ F i2 



\y\ ai 



\uo(y)\dy. 



Let us introduce the weight function w(y) — (— Pmr~ ) and the weighted L p norm 

,i/p 



IMU, 



(J H , Kw)IMv)di/) /p , if 1 <?<+«>, 

[ ess sup yeR3 \v(y)\, if p = +oo . 

L 2 conservation and (16.31) yield the estimates 



«(•,*) 




U 




«(',*) 


y/w 


LI, 


\/w 


7 

LI, 


VwJ 



^c 



l + t Q1 



ti 



(/o 



if,. 



Then, letting, for all p > 2, 



9 P = 1- 2 -, P ' = ^, 
P P~ 1 



so that 



9*6(0,1), 

■pola 



1 1 



p' 2 

the Riesz-Thorin interpolation theorem yields 



9 1 1 
1' p-~ 



<C"*> 



1+i" 1 



IIQ 



i.e. 



Hy,t)\ p C-^ L ) 2 ~ P dy) 1/P 



^C 



ti 



\u (y)f C-^) 2 - p ' dy 



w 



Hence, inequality (|1.29|) in Theorem 11.111 follows . Therefore, in order to prove the 
theorem, it is sufficient to prove estimates (I6.1[) and (|6.2j) . 

It is well known that the link between plane waves and a combination of zonal 
functions is given by the Jacobi- Anger expansion, combined with the addition the- 
orem for spherical harmonics (see for example |51j . [39j and the references therein). 
For N = 3, we get that 



(6.4) 



OO 

^E^(Miyi)4/ N (2//M)- 

£=0 



We need to estimate the kernel K in (|1.15|) which, as observed in (|1.27|) . can be 
written as 

K(x,v)=s(\x\\y\,ft,ft) 
where 

OO 

1=0 
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\ - y/(£ + l/2) 2 + a. We split the sum into two 



with be = y/(e + l/2) 2 + a, o 
terms 

(6.5) s(r, e,e') = J2 rbe i-^ ^ z o e) W) + E rbe i-»* Wo* &) 

t=o e=e 

= s 1 {r,e,9') + s 2 {r,e,e'), 

with -^o^O such that ag > for any £ < £$ and ag < for any € ^ £o (^l is 
meant to be zero if Iq = 0) . Our goal is to show that the singularities in the Jacobi- 
Anger expansion are described by the first finite sum Si at the right-hand side of 
(|6.5p while the second term S2 at the right-hand side is uniformly bounded. Such 
boundedness for S2 follows from the arguments below. We have that 

CO 

(6.6) S2= y £r bt j-a,{r)Z^{9 f ) 

t=to 

= £ i-^hWf) {9') + f;«-^i_«(r) - C^je(r))Z^ {d >) 



t=(\, 



= l 2 



(27r )-t e -^' - 2 r^wzf (*') 



00 

The first term at the right hand side of (I6.6[) is clearly bounded, since it is the 
difference between a plane wave and the first (£0 — 1) terms of its Jacobi- Anger 
expansion. 

We first notice that the second term at the right hand side of (|6.6[) is bounded 
for r ^ S if S > is sufficiently small. Indeed from the estimates 

(6.7) |J„(t)| < -J—^ty e t2 /\ for all „ > 0, t> 0, 

|Z^ (0') I < Z< £) (0) = ^ti, for all * > 0, 0, 0' G § 2 , 



see for example 39 , it follows that, if r ^ 5 

00 00 

Ez-S- a »4%) <E 

e=e e=£ 



2£+l (|) bg r2/4 



4tiT(& £ + 1) r \ 



C 



E 



V24.^ o r ( 6, + i) 



I)© 



-at e 



*74 



V2 4tt 



(I) "E 



2£+l 

r(6, + 1) 



sC CY _Q <o 



for some constant C > dependent on (5 and £q but independent of ?-, 0, 0'. 
Next, for r > 6, we write 



£(rS_a,(r)-rC<+*>i<(r))zJV) 



=4> 



i/ ; >(-«|N«^-. 



2 



w, 
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where we have used the following representation for Bessel functions 

:i f n( z -±ydz 



Mr) = jr-7 / e 



2wi J 1 z u+l ' 

with 7 being the positive oriented contour represented in Figure[T](see [351 5.10.7]). 
We have also exchanged sum and integral, which is allowed for any r,9,9' , as we 
will see below. 

For convenience, we split the integral along 7 into the integrals I\, along the 
circumference of radius 1 (to be denoted as Ti), and the integral I2, along the lines 
running between z = —00 and z — —1 (to be denoted as T%): 



Notice that, by analyticity of the integrand outside z = B.~ + i0 ± , we can write 




where T\ is the circumference of radius 1 + e around the origin and T| runs along 
(—00, — 1 — e) + iO . Notice that, for the integral along T\ U T| since \z\ > 1, one 
has absolute convergence for any given r, 9, 9' and hence the exchange of integral 
and sum performed in formula (|6.8[) is allowed by Fubini's Theorem. 
We start estimating the integral along Ti. Taking into account that 



we have that 

a 2 (log( ? z)) 2 , 0(1) 



(6.9) (izY t+ *>- be - I = logfc, 

y ' L v ' J 2^+1 sv ; 2 (2^+l) 2 e 3 

= Ji,i(z,£) + Ji,2(z,l) + Ji,z{z,£) 

as £ — > +00 uniformly with respect to z € Ti. Since z~ hl and z^ + ^' have a branch- 
cut at z € R~ , the function \og(iz) will also have a branch-cut at z 6 R _ , as well 
as the function (iz)^ that will appear below. From (J6.9P the contribution of the 
right hand side of ( 16.81) on V\ can be written 




(Ji,i(z,£) + Ji, 2 (z,£) + J h3 (z,£)) 



(iz) i+ ? 




J\,i + 1/1,2 + iTi 



1,3 



where every summand J7i.i, i = 1,2,3, corresponds to the integrand with the cor- 
responding Ji t i. Since on Ti we have that \z\ = 1 and then e 5 ^~~J = 1, from 
the estimate \zi (9') I < ^4^ we deduce that, if r > 6, 



in 



^ 21+1 



\J 1>3 \ < const r- 1 / 2 J2 -gT- < const (T 1 / 2 , 



fo 
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and hence |i7i,3| is bounded. Concerning 1/1,1, we notice that 



^ 



(«i-io.) ■"i..,(/:ix;^^w-*— ^^f;ft(«-«')«* 



£=o e=o 

1 a(iz)~2 log(iz) 1 o(— iz) 3 log(iz) 



where we have used the well-known identity (see for example [39 ) 

(6.11) 47r Ze (g,) = PA9-9'), 

v i 2£+l n h 

with Pi being the Legendre polynomial of index £, and the identity (see Formula 
22.9.12 in [2]) 

oo 1 

(6.12) Yp e (t)w e = - 

which is valid for \w\ < 1. Hence, identity (|6.10|) is valid for |z| > 1. Therefore 

(6.13) J,!- (Bounded terms) - -^ / e*(-*) to g(")H*)* Ji, 

V ^ V ; 87r 2 ir5 Vrj v / z 2 + 2 i z(0-0')-l z 

where the first term at the right-hand side of (|6.13p represents a finite sum of terms, 
that are needed to complete the series (|6.10p from i = to £ = £q — 1, and which 
are uniformly bounded. Since |eS( z-1 / z ) | = 1 for all r > and z £ Ti, if 1 — (9- Q') 2 
does not approach zero, i.e. if 9 ■ 9' stays far away from ±1, the second term at the 
right hand side of (j6.13[) is uniformly bounded with respect to e — > + , r > 5, and 
1 — (9 ■ 9') 2 > S, due to the integrability of the two square root singularities of the 
integrand at z± = -i{9 ■ 9') ± y/1- [9 ■ 9') 2 . 

When (9-9') = qpl, the two square root singularities at z± collapse into a stronger 
singularity at z = ±i. Let us discuss e.g. the case (9 ■ 9') = —1 (the case (9 ■ 9') = 1 
can be treated similarly); then 

(6.14) lim / e iH)Hi&M* 
e->0+ J r e z — i z 

2- ir r^Tr f Z( z -l)l0g(iz) . idz 

-K 2 ie lr + PV I e^ *> . (~iz) 2 — • 



p 6 — I 6 

Equation f|6.14[) is simply the Plemelj-Sokhotskyi formula (see for instance [T]) for 
the limit of Cauchy integrals when approaching a singular point. The hrst term 
at the right hand side of f|6.14[) is clearly bounded. The second term at the right 
hand side of (|6.14[) is a singular integral of the function e^( z ~~) og ( lz '(~ lz > — which 
is differentiable for z = e %e with 9 in the neighborhood of ? (remind that the 
discontinuity of the argument of z is along the negative real line) . 

Hence, since the principal value of a Cauchy integral of a differentiable function is 
bounded (cf. p]), we conclude the boundedness of J\,\ for any r > S. The fact that 
the principal value integral is bounded for any r does not exclude the possibility 
of its diverging as r — > oo. In order to exclude this possibility, we consider a 
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neighborhood in T\ oi z = i: 

r*° = [z = e*(* + *), | s |<s «l} 
and we integrate there for any r»l having into account that 

(iT + s)e J 2 7T 



i(e' ls — 1) s 

Hence, 



0(1) for s ~ 0. 



( " ~* s '"'~ '■- "- t - ircossl n ' 0(1)) ds 



i(e ts — 1) 7 s 

-so 
so 



n- ds + PV / O(l)e lrcoss ds. 



-so 



so 



Since PF J 7r- ds = 0, it follows that 



-so 



7 i{e ls - 1) 



-so 

"ISO 



Hence, the integral along T^° is uniformly bounded for any r. The integral over 
r^I^ is also uniformly bounded since e^V 2- */ = 1. Hence, the principal value 
integral over Fi is uniformly bounded. 

If one considers the two singularities z± sufficiently close, then, similarly to 
(|6.14p . the integral over fi can be written as 



(6.15) 

'Ti J Arc(z-,z + ) vTi\Arc(z„,z + ) 

where Arc{z-,z+) is the small arc of Ti between z+ and z~. The second integral 
at the right hand side of (16.151) can be easily estimated just like the principal value 
above and yields the same estimates uniformly in 8,6'. The first term at the right 
hand side of (|6.15[) is, after writing 

z ± = e<t±^) = _f(0 . Qi) ± y/l-(0.e') 2 , 

the integral of 

e lrcoss (s + Tr)e l i 



l/v 



( s +f) - A*-*)} (e^+i) - e < ¥+ f ) 



- 1 , if s < — s 

:(1 + 0(S-S) + 0(S + S)) X { +1 lis>S 



^s + s\\s-s\ ,..,._ 

v ■ —i if - s < s < s 

7re ircoss I I \ / 1 



y/\s + s\\s-s\ \y/\ s -s\J \V\s+J\ 
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where \v(s)\ = 1. As a consequence, we can estimate 

7re 

I Arc(z^ ,2+) 



7TP 



ds + O(l) 



s — s^ 

1 ircos(si) 



-rfi + 0(l) 



Li VT^ 

^ const, 

uniformly with respect to r and s. Therefore, we conclude that the integral on Fi 
is uniformly bounded both in S and r. 

Finally, the term Ji,2 in (|6.9|) . inserted at the right hand side of (|6.8|) . produces 



2nir 
where the series 



1_ r , (z _, A a\\og{iz)f (^ Zf(ff 



3 (w)- 



2; 



?W ^n 



7 W /-flM 



e=o 



(2£+l) 2 



9',z)dz 



is the primitive in z of the series 
g(6,0',z) = 



00 7W cflM 



e=o 



U+\ 



Thus, using (|6. 11|) and (J6.12I) . we conclude that F(0,8',z) is the primitive of 



9(9,9', z) 



(iz)- 



in ^1 + 2(9- 0')l-?r 



and hence, since g(9,9',z) presents a square root singularity if 9 ■ 9' 7^ —1 or a 
1/(2; — i) singularity at z = i \i 9 ■ 9' = — 1, we conclude that F(8,8', z) presents 
at most a log-type singularity, which is integrable, and consequently the integral 
yields a uniformly bounded contribution J\p,. Therefore, we conclude that I\ is 
uniformly bounded. 

We continue estimating I2, 



h 



1 



2nir2 jp 



ei(*-i)(£[(zz 



>e+h-bt 



1 



^TO 



(«) 



t+h 



dz 
z 



Introducing the changes of variables, z = e ±7 "e*, exchanging sum and integral 
(arguing as above) and rearranging terms, we rewrite it in the form 



1 00 



1 



where 



and 



2nir2 



Mr)=- 2f f bl { i!* ) I < ■-"• 



2ixir 2 



r(^2,i + ^2,2) 



ht/ e -6/* _g-(^+|)* 



)dt 



B<(r) = -2t / ■ 



-rsinht -(£+i)tf SJ1) "''' 



r/f. 
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We estimate Ag by using again that | Z t 



(f) 



< 21+1 
** 4,-k ■ 



\Jl, 



J2z^(e')M 



=*0 

oo 



^C]T(2^+1) / 



-r sinh t 



e _ M _ e -(^+|)t 



rft 



^Cj](2^+l) 



&£ £ + 



£C. 



In order to estimate Bg, notice that 

'sin7r&£ sin7r(^ + |) 



i'+* 



1 



l ' 2 

(-1) £ 1 / an 



sinnbg — sm7r(i + — ) 



sin7r6 



l ' 2 



(-iY 



(ITT 



i£+i2\2£+lj »*+* ^2(2^+1) 8(2^+l) 2 
i^i an bi 2 



o(r 3 ) 



an 



21+12 8 \(2^ + l) 

By using formula f|6.12[) . it readily follows 



+ 0(£-*). 



\J2,2 



J^Z^(9')B e (r) 



< 



e' rsinht e-^x 



K io (r,6,6',t)-il 



1 



-E 



° 5o 2 7r 2 t / -Je- A zi /) "»' 



4 ^/i _ j e -*(0 • 0') - e- 2 * 
) f 1 



16 



(2^+ir 



o(r d ) 



where J<Q accounts for the terms that need to be added in order to use (|6.12[) and 
which is uniformly bounded. 

Since ^/|1 - ie~ t (9 ■ 9') - e~ 2 *| ^ co i " st ^* for some const > and, using again 

that \zf {9') K *&, 



E * Vft ^ (') 



i 



(2^+ir 



■o(r 



°° e~ ft 



^C^ — ^2Clog|i| 



for some constat C and any 9, 9' , we conclude the existence of other constants 
C',C" such that 

1 + Vi 



s~V poo 



lht. 



Vt 



log |t| 



<<7". 



Hence the uniform boundedness of ^2,2 follows. We conclude then 
(6.16) sup \S 2 (r,9,e')\ < +00. 

If a\ = ao < 0, then £q — 0. Hence Si(r, 0, 0') = and (|6.1|) is proved. 
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If o.\ — ao > 0, then £q > 0. From (|6.7|) and the fact that ae ^ clq for all £ G N, we 
deduce that 

(6.17) |Si(r,0,0')| < const r _ao = const r"" 1 for all r < 1, 6», <9' g S* -1 . 
On the other hand, from (J5.4I) and (|5.5[) we easily deduce that 

(6.18) \Si (r, 0, 0') |< const for all r ^ 1, 0,0' e S^ -1 . 

Estimate (JO) follows from ([6"3]l . (|63g|l . ([637) . and ([B38"]) . D 




Figure 1. Integration oriented domain 7. 



Appendix 

The asymptotic behavior of eigenvalues /i&(A,a) as fc — > +00 is described by 
Weyl's law, which is recalled in the theorem below. We refer to [41] [44] for a proof. 

Theorem A.l (Weyl's law). For a € L°°(S A '- 1 ,]R) and A € C^S^ -1 ,^), Zet 
{/ifc(A, a)}fc^i oe ifte eigenvalues of the operator La. a = ( — * Vgjv-i + A) + a(0). 
Then 

(A.19) A*fc(A,a) =C(iV,A,a)fc 2/(Ar ~ 1) (l + o(l)) as fc -> +00, 

/or some positive constant C(N, A, a) depending only on N, A, and a. 

The following lemma provides an estimate of the L°°-norm of eigenfunctions of 
the operator La, a in terms of the corresponding eigenvalues. 

Lemma A.2. For a e L°°(S Ar - 1 ,R), Ae C 1 {S N - 1 ,R N ), andfceN\{0}, let tjj k 
be a L? -normalized eigenfunction of the Schrodinger operator La, a on the sphere 
associated to the k-th eigenvalue /Xfc(A,a), i.e. satisfying (I1.12[) . Then, there exists 
a constant C depending only on N , A, and a such that 

|^(*)|<chk| L(JV - 1)/2J , 

where \_-\ denotes the floor function, i.e. \x\ := max{j eZ: j ^ x}. 

Proof. Using classical elliptic regularity theory and bootstrap methods, we can 
easily prove that for any j e N there exists a constant C(N,A,a,j), depending 
only on j, A, a, and N but independent of k, such that, for large k, 



Uk\ 



< 



W '(JV-i)-2U-i) (S™- 1 ) 



C7(iV,A,a,j)(^(A,a)) 3 
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Choosing j = I Z 1 \ , by Sobolev's inclusions we deduce that 

for any < a < 1 — N ~ l + I jv 9 ~ 1 j , thus implying the required estimate. D 
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